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A higher-order bending theory is derived for laminated composite and sandwich beams thus extending the recent 
{ 1,2 1 -order theory to include third-order axial effects without introducing additional kinematic variables. The 
present theory is of order (3,2) and includes both transverse shear and transverse normal deformations. A closed- 
form solution to the cylindrical bending problem is derived and compared with the corresponding exact elasticity 
solution. The numerical comparisons are focused on the most challenging material systems and beam aspect ratios 
which include moderate-to-thick unsymmetric composite and sandwich laminates. Advantages and limitations of 
the theory are discussed. © 1998 Published by Elsevier Science Ltd. All rights reserved 


INTRODUCTION 

Higher performance and lower cost requirements for the 
next generation of aerospace vehicles often necessitate the 
use of advanced polymer-matrix composite materials. 
Composite materials can be tailored into highly efficient 
structures that combine high stiffness and strength, light 
weight, and improved fatigue and thermal performance. 
From the design perspective, accurate strain and stress 
predictions are required to avoid higher factors of safety that 
inevitably lead to an over design, reduced performance, and 
higher cost. 

The structural modeling of composite and sandwich 
laminates with the use of approximate beam, plate, and shell 
theories is known to be the most efficient. Many significant 
developments in this area can be found, for example, in 
review papers by Reissner 1 , Reddy 2 , and Noor and Burton 3 . 
Numerous finite elements used in commercial and research 
codes have also been developed for composite structures. 
The most commonly used finite element models are those 
based on the first-order shear-deformation theory. The 
following brief discussion reviews the most pertinent 
aspects of composite beam theories, and also makes 
comparisons to similar plate and shell theories. 

The classical Bernoulli-Euler beam theory, neglecting 
transverse shear and transverse normal deformations, is 
appropriate for thin, homogeneous beams and is known to 
be inadequate for composite and relatively thick beams. 
Timoshenko beam theory includes transverse shear defor- 
mation and provides more accurate response predictions for 
thin and moderately thick homogeneous beams. Reissner 4 
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stress-based and Mindlin 5 displacement-based first-order 
shear-deformation plate theories were originally developed 
for the analysis of homogeneous elastic plates. Many 
subsequent shear-deformation theories, utilizing the 
displacement-based approximation approach, focused on 
the analysis of laminated composites, e.g. refer to Stavsky 6 , 
Yang et al . 7 , Whitney and Pagano 8 . Such theories, 
commonly referred to as single-layer theories, treat a 
laminate as an equivalent single layer, with the displace- 
ment assumption representing a weighted-average distri- 
bution through the thickness. 

Reddy and Liu 9 formulated a layer-wise theory which 
assumes piece- wise smooth displacement components 
through the thickness, i.e., while the displacement function 
is continuous through the thickness, the slope of the function 
at the ply interfaces may not be continuous. This type of 
theory produces a large number of unknowns and is 
computationally expensive, especially when a laminate 
consists of many layers which is usually the case in load- 
carrying structures. 

Higher-order theories, which account for transverse shear 
and transverse normal stresses, generally provide a reason- 
able compromise between accuracy and simplicity; how- 
ever, they are usually associated with higher-order boundary 
conditions that are difficult to interpret in practical 
engineering applications, e.g., refer to Essenburg 10 , 
Whitney and Sun 11 , Lo et al . 12 , Reddy 13 , and Phan and 
Reddy 14 . 

Recently. Tessler and coworkers 1 '' 19 developed a 
higher-order theory for application to laminate composite 
beam, plate, and shell analyses. The theory maintains the 
simplicity and computational advantages of the first-order 
shear-deformation theory. It accounts for transverse shear 
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and transverse normal deformations by assuming a special 
form of the {1,2} -order displacement assumption (The 
notation { mji } implies that the axial displacement is 
expanded with a polynomial of degree, m, whereas the 
transverse displacement may be of a different degree, n.) 
Additionally, the average shear strains are assumed to be 
parabolic, thus satisfying zero shear tractions on the top/ 
bottom surfaces; and an average transverse normal strain is 
assumed in the form of a cubic polynomial satisfying one of 
the equilibrium equations of elasticity theory exactly. The 
approach requires that the transverse strains need only be 
least-squares compatible, through the laminate thickness, 
with the strains derived from strain-displacement relations. 
The resulting thickness distributions for the transverse 
stresses and strains produce adequate correlation with 
results given by elasticity theory, an improvement over 
previous higher-order theories. Tessler 16 improved the 
theory further for application to composites by introducing 
an independent polynomial assumption for the transverse 
normal stress to replace the cubic transverse normal strain 
assumption. The improved theory results in a more accurate 
representation of transverse normal stresses and strains, and 
is further substantiated by solutions given by Schleicher 20 . 
The {1,2} theory retains the simplicity of the first-order 
shear deformation theory in so far as the engineering 
boundary conditions are concerned. Furthermore, the 
theory gives rise to finite element formulations that are 
fully compatible with the first-order shear deformation 
elements. 

Application of the {1,2} theory generally results in 
excellent predictions for thin and moderately thick homo- 
geneous and laminated composites. Nevertheless, the theory 
has some limitations, particularly with respect to the 
modeling of relatively thick sandwich laminates. This is 
because in such laminates the distribution of the inplane 
displacement and strain can be highly non-linear. In thick 
laminates, this generally results in underestimation of the 
axial stress, typically the largest stress that governs the 
design of the structure. Another deficiency, which is only 
manifested in sandwich laminates, is the violation of the 
traction conditions on the top and bottom surfaces 
associated with the transverse normal stress. 


In a recent NASA publication, Cook 21 explored a {3,2}- 
order >eam theory which expands upon Tessler s {1,2} 
theory by including cubic axial effects. A special hier- 
archic; 1 form for the axial displacement is developed such 
that th? theory employs the same five kinematic variables 
as its {l,2}-order counterpart, without introducing any 
additional kinematic variables. The hierarchical form of the 
displacement field ensures the exact fulfillment of traction- 
free snear stress boundary conditions and permits a 
straightforward reduction to several lower-order beam 
theories. As in Tessler 16 , in addition to the assumed 
displacements, an independent polynomial expansion is 
ernplo; r ed for the transverse normal stress. The concepts of 
transvt rse shear and transverse normal correction factors 
are effectively incorporated using strain energy and traction 
equilibrium considerations. The theory enables more 
accura e predictions for the axial, transverse shear, and 
transvt rse normal stresses and strains, particularly for thick 
laminated composite and sandwich beams. Accurate piece- 
wise smooth transverse shear stresses are determined by 
integrating two-dimensional equilibrium equation of elas- 
ticity theory. Cook 21 also developed a straightforward 
correction procedure that improves the accuracy of this 
approach for unsymmetric and sandwich laminates. 

In this paper, the theoretical foundation and predictive 
characteristics of the {3,2 (-order theory are closely 
examined. The theory, which begins with an assumed 
{3,2 (-order displacement field and assumed cubic trans- 
verse r ormal stress, employs the virtual work principle from 
which the beam equilibrium equations and associated 
boundary conditions are derived. These field equations are 
solved in closed form for the problem of cylindrical bending 
of laminated composite and sandwich beams. Appropriate 
transverse shear and transverse normal correction factors are 
employed. Numerical results are presented for moderately 
thick and truly thick beams, and comparisons are made to the 
{ 1,2} -order theory and three-dimensional elasticity solutions. 

{3,2}- ORDER BEAM THEORY 

Consic er a straight, linearly elastic beam laminated with N 
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orthotropic plies subject to the loading shown in Figure 1. 
The beam has a span L and a rectangular cross-section 
thickness of 2 h and width h. The orthotropic plies are 
stacked from the bottom (z = —h) such that the material 
properties, in general, are functions of the z coordinate. The 
tractions q + and q~ are applied normal to the top and bottom 
faces of the beam. and T lL (i = jc,z) are tractions 
prescribed at the ends of the beam. 

Displacement assumptions 

From the viewpoint of exact elasticity theory, the 
displacement components are piece-wise smooth and, in 
thick laminated composite and sandwich beams, they are 
non-linear through the thickness. This contrasts with the 
predominantly linear displacement distributions for thin 
beams. To represent both linear and non-linear deformation 
effects within the realm of a relatively simple, single-layer 
structural theory, the axial and transverse displacement 
components u x and u z are assumed to vary through the 
thickness as the cubic and quadratic polynomials 

3 

u x (x, z) = X M, Wf, u.(x, z) = w(x) + H'l (.r)f 

<t?, ( 1 ) 

+ vv 2W(r~ + o 


displacement components of the form: 

«,(*, -) = U + h& - im 1 - 1 )/m Lv - hS{ - i)Y, 
u.(x, z) - w + fvv, + (f 2 - ±)w 2 , (4) 

where 

y = + w.tW] + vv’ 2 , V (.v) 

Note that the resulting displacements in eqn (4) are in terms 
of the same five kinematic variables as in Tessler s {1,2}- 
order theory, and the quadratic u z is the same for both 
theories. The first three variables, m(x), 0(a), and w(x), are 
the Reissner weighted-average displacements, defined in 
eqn (2), whereas W]{x) and w 2 (x) represent the higher- 
order terms that account for the stretching of the beam 
through the thickness. The cubic u x has an hierarchical 
form such that if the higher-order terms nq v and y are 
eliminated, the displacement field is reduced to the { 1, 21- 
order theory with a linear axial displacement distribution. 

The displacements in eqn (4) represent the beam analog 
of the {3,2 (-order plate displacement approximations 
explored by Tessler 17 in the context of a hierarchical 
recovery of the {1,2} results using the {3,2} -order 
displacement, strain, and stress expansions. 

S tress -st ra in re l at i ons 


where f = zJh E [ — 1,1] is a dimensionless thickness 
coordinate such that f = 0 defines the midplane of the 
beam. The four w, coefficients in the axial displacement 
expression are yet to be defined, the w, coefficients in the 
transverse displacement represent the same kinematic vari- 
ables as those defined by Tessler 15 . The constant C is 
included in the transverse displacement equation to allow 
w(x) to represent a weighted-average transverse displace- 
ment yet to be defined. 

Three weighted-average kinematic variables are defined, 
as in Reissner 4 , such that 


u w=rhl-> u * (x ’ z)dz ' 6(x)= h 

3 f fl 

w(x) = — J f u z (x,z)(] - f 2 ) d Z 


u x (x 9 z)z d z. 


( 2 ) 


where u(x) is the midplane displacement along the x axis, 
0(a) is the rotation of the normal about the y axis, and w(x) is 
the weighted-average of the transverse displacement. The 
displacement field in eqn (1) is substituted into eqn (2) 
resulting in C — — V 5 , and the expressions for two w, 
coefficients in terms of the m(a) and 0(a) variables. The 
remaining two u, coefficients are determined by enforcing 
zero shear traction conditions at the top and bottom beam 
faces. Since, from Hooke’s law (eqn (5)), the shear stress is 
proportional to the shear strain, the shear strain at the top 
and bottom faces must also vanish: 


Ixz = 



(3) 


Enforcing the conditions in eqn (3) gives rise to the 


Either plane strain or plane stress constitutive relations 
can be developed for a laminated beam resulting in the 
stress-strain relations in the form: 
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The complete derivation of eqn (5) can be found, for 
example, in Cook" 1 . 


Strain -displacement relations 

In this beam theory, two distinct approaches are used to 
express the strains in terms of the kinematic variables. The 
axial e xx and transverse shear Y vr strains are determined 
directly from the strain-displacement relations of elasticity 
theory, where the y xz strain is further augmented with a 
shear correction factor. The strains derived in this manner 
will be represented by continuous and differentiable 
polynomial functions through the thickness whose distribu- 
tions are independent of the individual laminate properties; 
hence the superscript (&) will be dropped for these strains. 
Clearly, such strains can be regarded as some average 
representations of the "true” strains (i.e., those strains that 
satisfy the requisite equations of elasticity theory). The 
derivation of the transverse normal strain e[ k : \ however, 
begins with an average stress assumption for o :z , which is 
assumed to have a cubic polynomial distribution through the 
thickness. These strain developments are summarized as 
follows. 

The axial, average strain e xx is obtained from the linear 
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strain-displacement relations as 

C.V.T = «v,.v = + %\i)4* [ + + K H <t > 3 (6) 

where the strain measures, curvatures, and the thickness 
distribution functions <£, are defined as 

(e >o.£w) = (m. v ,Aw| > .„), 

(k<<). *tf) = (0. V u '. -v.v + 0.,) + m '2,-u)' (7) 

(4> u <t>z*<h) = (*r, i/6 - r 2 / 2 , /i(f/5 - r'/3)) 

The transverse shear average strain is obtained from the 
linear strain-displacement relations of elasticity and is 
augmented with a shear correction factor fc, i.e., 

yT ~ k y<: = *(«*.: + «-..«) = 

, , ( 8 ) 

(Tv^)' <Av.-) = (0 + v 5( 1 - ?')/4) 

The shear correction factor is introduced in eqn (8) in antici- 
pation that for certain material systems and lay-ups, a 
correction in the value of the transverse shear energy may 
be necessary; the shear correction factor provides a simple 
and effective mechanism for implementing such a 
correction. The motivation for circumventing the determi- 
nation of the transverse normal strain directly from the 
strain-displacement relations is as follows. The strain- 
displacement relations which employ the displacement 
assumptions eqn (4) give rise to a continuous through the 
thickness e zz strain which would represent only an average 
distribution of this strain through the thickness. This in turn 
would result in a <j :: which for laminated beams may exhibit 
discontinuity along ply interfaces. However, according to 
elasticity theory, a :: must be continuous through the thick- 
ness and el* 1 may be discontinuous along ply interfaces. The 
approach introduced by Tessler 16 enables the derivation of 
an improved el* 1 that will be discontinuous at ply interfaces. 
Importantly, the desired simplicity of the theory is retained. 
For mechanical loading, o :: is closely approximated by a 
cubic expansion through the thickness as 

3 

°::= X a -»r 
«=() 

in which the four a zn coefficients need to be determined. 
Two of the coefficients are found from the equilibrium 
equation of elasticity theory, i.e., 

t vc, x “f - = 0 (10) 

Since the transverse shear stress satisfies traction-free 
boundary conditions on the top and bottom surfaces of the 
beam, i.e. r A: (.v, ± h) — 0, the derivatives of the shear stress 
r XZ K at the top and bottom faces must vanish. To satisfy the 
equilibrium equation, the derivatives of the transverse 
normal stress must also vanish on the top and bottom 
surfaces: 

±A) = 0 (11) 

These exact equilibrium traction conditions reduce o :: to the 
form 

°:z = tf-0 + <7; 1 05* <t>5 = (f “ ( 12) 


The remaining two coefficients are found by forcing the e l* 1 
strain :o be least-squares compatible with the corrected 
averag ; strain derived from the strain-displacement relation 
(the notation el k J with the superscript ( k ) implies that the 
strain is piece-wise (ply-level) continuous): 

rh 

minimize J f (el* 1 — d z (13) 

where the corrected average strain is determined from the 
strain-displacement relations as 

u L - . — k~n)E z o + 2k z \ K z i)<t >\ , (c^o» %>) = (vv j //i, H r 2 fh ) ( 1 4) 

where k z{) and k z \ are the transverse normal correction fac- 
tors, ai d £, 0 and k^) denote the transverse strain measure and 
curvature, respectively. Obtaining el* 1 from the constitutive 
relations, eqn (5), results in 

eW = l/c£V:, - C\ k Js xx ) (15) 

Introducing eqns (14) and (15) into eqn (13), where the 
minimization is performed with respect to the undetermined 
coefficients, a# and o zU results in two algebraic equations 
from which these coefficients are readily determined. Eqn 
(15) is then simplified to yield the transverse normal strain 
of the form 


= e m 4>\ V + k :0 e ^)'p t 2 ) + e H^T + K.xO'I'T 

+ k zl k 3)'I / 5 ' + ' 


(16) 


where depend on the thickness coordinate, and the 
elastic stiffness coefficients, C^l - For their explicit form, 
refer to Cook 21 . 

In contrast to the linear distribution of in eqn ( 14), 
el* 1 ca l be discontinuous at the ply interfaces and is piece- 
wise cubic. As will be demonstrated by numerical 
comparisons with exact elasticity solutions, this form of 
el* 1 ensures superior through-the-thickness predictions and 
improves the overall beam response. 


Variat onal principle 

The principle of virtual work is employed to construct the 
beam equilibrium equations and associated boundary 
condit ons. Neglecting body forces, the virtual work 
principle can be stated as 

J'^(a** , Se rT + a^Sel* 1 + 7^67 xz ) d A dr 

— q ' 5u z (x, h) dv dy + q &u z (x, — h) dv dy 
4- I [r rf) «M,(0^) + 7\ <) «*i-(0,z)]dA 

' A 

- [T xL bu x (L, z) + T :L 6u z (L, z)] d4 = 0 (17) 

I A 

where 6 is the variational operator, A is the cross-sectional 
area o ’ the beam, and S + and S~ denote the top and bottom 
surfaces of the beam, which, respectively, are subject to the 
normal pressure loads q + and q ~ . The first term in eqn (17) 
is the olume integral representing the virtual work done by 
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the stresses. The surface integrals denote the virtual work 
done by the external surface tractions. 

Introducing the beam displacement assumptions and 
strain-displacement relations into eqn (17), then integrating 
over the beam cross-section and performing integration by 
parts results in the one-dimensional form of the virtual work 
principle 

£ [N x x du+{Q,-M x ,,-- 4 MH.M e 

+ ( 4 — Q xx — q \ )Sw + (NJh + hN H xx — < 72)^1 
+ ( MJh 2 + M H xx — )5h ! 2 ] dLv 

+ X (l-2a/L)([^ ra -/V f (a)]6a(«) 

c* = 0,/. 

+ [M xa - M x (a)}86(a) + [M Ux - hN H (a)]&w Ux (a) 

+ [M 2a ~ M H (a)]5{ f [0(a) + w’, A (a)] + w’ 2 . ( (a)} 

+ [Q.x* - Q x (u)}&w(a) + [Ola - hN Hx (a)5w { (a) 

+ [02a - M Hx (a)\ bw 2 (a)) =0 (18) 

where the beam reactive and prescribed (superscribed with a 
bar) stress resultants are defined as 

{N x ,N,,N h )= \ a (a ( *> + a. : <\ k*o zz if\ 0^2 

+ <7~\^3 ^ d/4, 

Q,= [/*’?>,, dz, {M x ,M.,M h ) 

= {4 (<&’*. + *:l Vzz'I'V’ dxxfa 

+ a z:'Pb ) ) dA - 

(q u q 2 )=b(q + -q ,q + +q~), 

A\a - T m dz, M xa = \j xa z dz, 

M\ a = h ^T xa 4> 2 dz, M Ia = ^T ra <h dz, 

Qxa = J ^ T 7a dz, 0 1 a = T- a f dz, 

02a = T ra (f 2 - I)dz, (a = 0, L) ( 1 9) 

Equilibrium equations and boundary conditions 

The equilibrium equations and boundary conditions are 
obtained from the principle of virtual work, eqn (18). The 
expressions associated with the arbitrary kinematic varia- 
tions must vanish independently, resulting in the following 
equilibrium equations: 

AC., = 0, NJh + hN H xx — q 2 — 0, \M h xx — Q xx — q { =0, 
0, - M,., ~ v = 0. MJh 2 + Mff.xx - % =0. (20) 

The remaining terms in eqn (18) must also vanish 


independently, thus giving rise to the boundary conditions 
for the theory. Evidently, either tractions or displacements 
can be prescribed at .v = 0 and L, such that 

N X0l — N x (a) or 8u(a) = 0, M xa = M x (a) or 80(a) = 0, 

M] a = hN H (a) or 8w\ v (a) = 0, A/ 2a = A ^h(°0 or 67 (a) = 0, 
Qxa = QxM or 5 w(a) = 0, Q\ a = hN H x (a) or 6 vv,(a) = 0 , 
Q 2a = M H Ja) or 6 vv 2 (a) = 0, (a = 0 ,L). (21) 

Beam constitutive relations 

Introducing the strains eqns ( 6 ), ( 8 ) and (16) into the 
stress-strain relations eqn (5), and integrating the stress 
resultants eqn (19), yields the beam constitutive relations of 
the form 

'N x ' 

N : 

N h 

< M, > 

S4 : 

M h 

, Q, > 


^11 

C*0^12 

^13 

Bn 

Cl *12 

B 13 
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k:\A \2 

k j/t 2 2 


k*B 2 , 

k~[)k. | ^22 

C <)^23 

0 


Co A 2 3 

A 33 

Bn 

C1C2 

B» 

0 

Bn 

C0C1 

Bn 

On 

*-,0 12 

On 

0 

C1C2 


C 1 By 2 

k~\ D \2 

k:,D 22 

Ci O2} 

0 

Bn 

C-0^23 

Bn 

On 

C 1 o>i 

On 

0 

0 

0 

0 

0 

0 

0 

k 2 G 


K .-o 

e H 
*>i) 

*z<> 

k H 
7vr(l 

where the A yj , B ip D ,, and G coefficients represent the mem- 
brane, membrane-bending coupling, bending, and shear 
rigidities. For their explicit form, refer to Cook 21 . 

Equilibrium equations in terms of displacements 

To facilitate a closed-form solution to the equilibrium eqn 
(20), subject to the appropriate boundary conditions eqn 
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( 21 ), it is convenient to express eqn ( 20 ) in terms of the five 
kinematic variables of the theory. First, the strain measures 
and curvatures can be expressed in terms of the kinematic 
variables in matrix form as 


e A<) 


~ it 
i lx 

0 

0 

0 

0 ■ 


e c0 


0 

1 

/» 

0 

0 

0 


/ \ 
u 

£ h 


0 

dX~ 

0 

0 

0 


VV] 

< *xi) 

> = 

0 

0 

0 

il 

fix 

0 

< 

vv 

4-0 


0 

0 

0 

0 

1 

fir 


e 

k H 


0 

0 

5 <r 
4 H\- 

5.JL 
4 fix 

'r 

iKx~ 


< W 2 J 

< Tvco j 


_ 0 

0 

<Kx 

1 

0 _ 



Substituting eqns (22) and (23) into eqn (20), the equili- 
brium equations in terms of the kinematic variables take 
the form 


1 3 kxxx + Jh ' 1 , w + W | 4 v 4 D l3 6 xxx 

^ W 2,xx + All [K^.va.v 4 h Ua) + w 2Ax\ 

” M , A 3 “ ^ 1 +k z I hB^W I , c.t 

/ n a I ^" 1 ^ 0^22 

— K-|Dp0 t H — --j W 2 

n 

- - k z] D 2 3 [ f (0 A 4 VV o) 4- w 2 ,xv] | “ = 0 (28) 


where v Ax = u* uu . 

Similarly, the boundary conditions eqn (21) can also be 
readily expressed in terms of the five kinematic variables if 
necessary. Eqns (24)— (28), subject to the boundary con- 
ditions eqn ( 21 ), can then be solved simultaneously to 
determ ne the five kinematic variables and subsequent 
displacement, strain, and stress distributions in the beam. 


Reduction to lower-order theories 



k^)A j 2 




4 

h U| - 

v 4 hA^w Uxxx 4 B u u m vv 



k - j B j 


5, x 


4 

c w '- 

v 4- B\i 

^(^. VA 4 W xxx ) 4 VV 2( ya y 

= 0 





(24) 


The lierarchical displacement approximation of the (3,2 } 
order theory permits a straightforward reduction to several 
lower-order theories. By eliminating the higher-order 
displacement terms h’i v (a) and -y(jc) from eqn (4), the 
displacement field reduces to the { 1 , 2 } form given by 
Tessler 1 


k ~()4 ■> ^ 

^ 1 3 11 . AAA 4- V\ | v v 4 hA 3 3 H | 4 v 

h 

4 ^ 31 ^. AAA 

k B -> 

+ ^2 VV 2 . VA 4- #33 [ 4 (® t yaa + vv . Ax 

) + 42,4a] | 

+ ^ ^ k A} A 1 2 u v 4 ^ w j 4 h vv j v r 4 ^-0^2 1 0 x 

k^4)k~ I B~n , _ r « 

+ U 2 + ^ - 0^23 [ 4 (0 V 4- w 

vv) 4* W 2 vv ] | 

-r /2 =0 

(25) 

S f D , ^: 0^23 , , D 

fi !3« AAA 4 H’ L u 4 hB^w { 

.4v + ^ 13 ®, AAV 

4- J t 2 U 2 . A V 4 D 33 [i(^.AAA + VV . 4v ) 4 VV 2 i 4y] | 


- k 2 G(6 x 4 vv vv ) -tj\= 0 (26) 

x.x + “ ^ “ * h’ i , , v 4- hB$ | w i AVV 4- D\ 1 0 vv 

+ '^2 1 VV 2,a + ^13 [ 1 ( 0 , va + 4aav) + vv 2,aaa] 

+ V y 4- — — W| , y 4- hBy$ w | vxv 4- D \ 3 Q xx 

+ vv 2, A + All + H’ y.ya) 4- H ? 2,xvx] | 

— k 2 G(6 + vv v ) =0 (27) 


u x (x\ z) = «4 «,(x, z) = vv 4 fw, 4 (f 2 - vv 2 (29) 

Consec uently, the higher-order strain and curvature terms 
e H and k h are eliminated from the theory. This results in the 
simplif cation of the equilibrium equations, boundary con- 
ditions and stress resultants, respectively, such that all of 
the terms with a subscript H are eliminated. 

The {1,2} displacement theory can further be reduced to 
Timoshenko theory by neglecting the Poisson effect (i.e., by 
setting >13 = 0 ), thus ignoring the coupling between the 
axial aid transverse stretching of the beam. Furthermore, 
the we ghting function associated with the computation of 
the trar sverse shear stiffness, which is parabolic, needs to be 
set to unity to simulate the constant shear distribution 
accord: ng to Timoshenko theory. While this yields the 
Timoshenko theory equilibrium equations, the boundary 
conditions for both {1,2 (-order and Timoshenko theories 
are the same. The results of Timoshenko theory can further 
be redi ced to those of classical beam theory by setting the 
transverse shear rigidity to be infinite, i.e., G = 

CYLIb DR1CAL BENDING PROBLEM 

The problem of cylindrical bending is considered for the 
beam n a state of plane-strain. The beam is simply- 
suppored at the ends x = 0 and x = L and is subjected to a 
transve rse load in the form of a half-sine wave applied at the 
top sur ’ace, i.e., 

</ + (a) = q„ sin( irx/L), q ~ {x) = 0 (30) 

where \ t) is the amplitude of the loading. 
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A closed-form solution is derived by first assuming 
appropriate trigonometric distributions of the kinematic 
variables 

u — U cos (nxIL), 6 = ( w ) cos(7ta7L), 
w = W sin(7 rx/L), vv j — sin(7r x/L), w 2 = W 2 sin(7 rx/L) 

(31) 

which satisfy the simply-supported end conditions exactly: 
At x = 0: N,( 0) = M x (0) = N H (0) = M H ( 0) 

= w(0) = vv, (0) = vv 2 (0) = 0 

At .v = L : N X {L) = M X (L) = N H {L) = M H (L) 

(32) 

= w(L) = w l (L) = w 2 (L) = 0 

Introducing eqns (30) and (31) into the equilibrium eqns 
(24)-(28) results in five algebraic equations in which the 
trigonometric functions are factored out, leaving only the 
amplitudes U , 0, W , W |, and W 2 as unknowns. Once 
the displacement amplitudes are determined, the kinematic 
variables are completely defined, giving rise to the strain 
measures and curvatures. The displacements, strains and 
stresses are then computed in a straightforward manner 
and are subsequently compared with the corresponding 
exact elasticity solutions, e.g., refer to Pagano" and 
Burton and Noor 23 . 

Since, in composite and sandwich laminates, the actual 
shear strain distribution is generally discontinuous at the ply 
interfaces and the shear stress is only piecewise continuous, 
the two-dimensional equilibrium equation of elasticity 
theory needs to be integrated to obtain an improved 
approximation for the transverse shear stress. This well- 
established procedure has been modified by Cook 21 to 
ensure accurate shear stress computations for unsymmetric 
and sandwich laminates — the type of laminates for which 
the integration approach results in rather inaccurate shear 
stresses. 


Numerical results 

Cook 21 assessed the {3,2} -order theory by examining a 
wide range of laminates and material systems. As expected, 


the best performance is achieved for homogeneous beams, 
where the displacement, strains and stresses, both due to the 
{3,2} and { 1,2} theories, correlate exceptionally well 
with exact elasticity solutions even for the thick beams 
with Lllh = 4. For homogeneous beams, all correction 
factors take on the value of unity (ft 2 = ft-o = ft-i = 1 . 0 ), i.e. 
no corrections are required. 

Presently, the numerical assessment is focused on the 
material systems and aspect ratios which expose the highest 
degree of modeling difficulty for the theory. In particular, 
the results for two types of moderately thick and thick 
composite beams (U2h = 10 and 4) are presented: (a) 
graphite/epoxy (GR/EP) unsymmetric laminated beams 
with a lay-up of [()4/90 4 /0 4 /90 4 ] T and (b) GR/EP, PVC- 
core symmetric sandwich beams with a lay-up of [O 4 / 9 O 2 /O 4 / 
9 O 2 /O 4 /PVC Core] s . The material properties and geometric 
data are summarized in Table /, and the transverse shear and 
transverse normal correction factors are given in Table 2. 
For details on the determination of the correction factors, 
the reader is referred to Cook 21 . 

Laminated composite beams 

In Figures 2 and 3, the displacement, strain, and 
stress through-thickness distributions for the moderately 
thick (U2h — 10) and thick (U2h = 4) unsymmetric, GR/EP 
laminated beams [ 0 4 / 90 4 /0 4 /90 4 ] T are shown. For compari- 
son purposes, the { l, 2 (-theory results are included for the 
thick case only where the differences in results are most 
pronounced. Due to the lack of symmetry in the lay-up, the 
midplane is in tension with respect to the e vv strain and is 
under a compressive e : . strain. The transverse displacement 
is non-linear through the thickness, and is within 0. 1 % of the 
exact solution for the U2h = 10 case and within 2% for the 
Li2h = 4 beam. The exact o :: stress is seen to be more 
complex through the thickness than its cubic approximation 
within the present (and {1,2}) theory. Nonetheless, the 
qualitative comparison is quite adequate. Also, the cubic 
distribution of the axial strain, e vv , is quite accurate, 
underestimating the maximum value only slightly for U 
2 h = 10. For Lf2h = 4, however, the results are significantly 
less accurate, with the present theory results being 
consistently superior to the {1,2} theory results. As 


Table 1 Material properties and lamina geometric data 
Graphite/Epoxy (GR/EP) 

Polyvinyl chloride (PVC) 

Ply thickness 
Beam width 

Sandwich core thickness 
Notation: L — Longitudinal direction, T = Transverse direction 


= 22.9 Msi 
E r = 1.39 Msi 
E= 15.08 ksi 


p LT = 0.32 
v n = 0.49 
v = 0.3 
0.00625 in 
l .0 in 

0.8 X Total thickness 


G ij = 0.86 Msi 
G n = 0.468 Msi 
G = 5.80 ksi 


Table 2 {3,2} and {1,2} theory correction factors 


Material system 

k 2 


k T <.) 


k :l 



{3,2} 

U,2} 

{3,2} 

{1,2} 

{3.2} 

{1,2} 

GR/EP laminate 

0.76187 

0.73262 

1.21668 

1.0 

1.01975 

1.0 

GR/EP-PVC sandwich 

0.30666 

0.37301 

1.24326 

1.24326 

1 .59569 

1 .76405 
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expected, the results are in excellent agreement with the 
exact solutions for the moderately thick beam and are 
somewhat less accurate for the thick beam. 

Sandwich beams 

Sandwich laminates present a unique challenge to any 
approximate bending theory owing to the drastic change in 
the material properties through the thickness. The face 


sheets of a sandwich are stiff while the core material is 
lightweight and, generally, is several orders of magnitude 
more compliant. 

Figures 4 and 5 show the displacement, strain and 
stress variations, through-thickness, for the moderately 
thick (iJ2h — 10) and thick (LI2h = 4) symmetric sandwich 
beams. Characteristically for a sandwich laminate, the axial 
stress c xx is carried by the stiff GR/EP face sheets whereas 
the transverse shear stress t x: is almost exclusively carried 





Figure 2 Unsym metric GR/EP laminate. L/2h = 10 
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by the PVC core. Note that for the moderately thick beam, 
the deflection is over estimated by about 2%. For the thick 
beam, the {3,2}- and { 1,2} -order theories over estimate the 
deflection by 12% and 37%, respectively. Such large 
discrepancies could have been avoided by way of correction 
factors appropriate for the thick regime. The axial 
displacement and strain have a pronounced zigzag distribu- 
tion through the thickness according to the exact solution. 
For these quantities, the cubic variations of the {3,2 (-order 







0 0.5 1 1.5 2 2.5 


theory predict the response adequately in the face sheets and 
at the midplane. Consequently, the stresses and strains on 
the top and bottom faces, where these quantities are usually 
the largest, are accurately predicted by the theory. Notice 
that the linear approximation for the axial displacement in 
the {1,2} theory underestimates the axial strain, resulting in 
a significant underestimation of the axial stress. The {3, 21- 
order theory captures <r vv at the top and bottom surfaces 
adequately, while the same stress for the {1,2} theory is 





Figure 3 Unsymmetric GR/EP laminate, U2h = 4 
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85% in error. The quality of these results suggests that the 
span to thickness ratio of U2h = 4 may constitute the 
practical limit for application of this theory to sandwich 
beams. 


CONCLUSIONS 

A { 3,2 (-order bending theory for laminated composite and 


sandw ch beams has been developed. The theory employs a 
hierarchical form of a third-order axial displacement and a 
quadratic transverse normal displacement, and possesses the 
same cinematic variables as the {l,2}-order theory. The 
assumed kinematic field results in an average parabolic 
shear strain such that zero shear-stress boundary conditions 
on the top and bottom beam surfaces are fulfilled exactly. 
An independent expansion for the transverse normal stress 
is also introduced, thus enabling accurate transverse normal 






T«/q 0 



°«/q 0 



Figure 4 GR/EP-PVC symmetric sandwich, U2h — 10 










strain and stress predictions. Appropriate transverse shear 
and transverse normal correction factors are used to adjust 
the shear and thickness-stretch response of the beam. A 
closed-form solution to the cylindrical bending of moder- 
ately thick and thick unsymmetric laminated composite and 
symmetric sandwich beams has been developed. The 
numerical results show that the {3,2} -order theory has 
some advantages over the { 1,2 {-order theory, particularly in 
predicting the axial response in thick sandwich laminates. 
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